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We study the kinetics of 2D Bose gas cooling provided Bose particles 
interact with 3D phonons. At low temperatures phonon emission is prohib- 
ited by the energy and the momentum conservation. We show that both 
particle-particle scattering and impurity scattering assist the Bose gas cool- 
ing. The temporal relaxation of temperature follows the law T ~ 1/ above 
the Berezinski-Kosterlitz-Thouless phase transition point and T ~ l/t after 
a Bose-Einstein 2D quasi-condensate develops. 

Exciton gas in GaAs bilayer quantum well represents a system where a 
2D Bose-Einstein quasi-condensation (BEqC) is possible at low temperatures. 
Experimental efforts ^ |2[ have been directed to produce such an exciton gas 
and to cool it down to a BEqC temperature. A short laser photo-illumination 
pulse excites electron and hole pairs. The so-called indirect exciton technique 
is used where a perpendicular electric field drags electrons and holes apart 
into two spatially separated layers. Then an electron and a hole bind them- 
selves into an indirect exciton particle. This experimental setup suppresses 
the electron-hole recombination giving rise to a relatively long exciton life 
time. During the initial photoexcitation pulse newly born excitons are hot 
and form a non-equilibrium state. In short time after the pulse ends the 
exciton gas reaches an equilibrium at some effective exciton temperature 
which is much higher than the lattice temperature of the cold GaAs crys- 
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tal. Frequent exciton-exciton collisions ensures the exciton temperature to 
be uniform across the bilayer. Exciton gas then start to cool down slowly 
due to an emission of phonons into a crystal away from the bilayer. This is 
the longest phase of the experiment limited only by a decay time of excitons 
due to the electron-hole recombination. In order to reach a BEqC point one 
needs both low temperature and a high density of excitons. Hence a fast 
cooling is essential. 

Recent calculation of energy losses in 2D ideal exciton system has pre- 
dicted an extremely slow cooling at low temperatures with the temporal law: 
T(t) ~ 1/ log(t), where t is the time 0. This fact is intimately related to the 
energy and momentum conservation which prohibits an emission of phonons 
by an exciton moving slower than the velocity of sound in GaAs crystal: c. 
Thus, the exciton gas cooling in effect stops when the exciton temperature 
falls below a characteristic blocking temperature: = mc^/2, where m is 
the mass of exciton, even if the crystal temperature is zero. 

This kinetic bottleneck problem becomes especially acute when the exci- 
ton gas is subjected to a strong perpendicular magnetic field that quenches 
the motion of exciton to the lowest Landau level and, thus, is helping to bind 
electrons and holes into exciton pairs. In this case the effective mass of an 
exciton is determined by the Coulomb interaction and can be much larger 
than either the electron or hole mass [Q. This results in a higher blocking 
temperature Tb and makes it difficult to reach low temperatures in the end. 

In this communication we supplement the analysis of exciton cooling of 
Ref. by an addition of exciton-exciton collisions and scattering on impu- 
rities. Both events substantially assist the phonon emission. We specialize 
to the case of exactly zero lattice temperature which allows us to neglect 
exciton-phonon scattering. Actually we are dealing with a general problem 
of 2D Rose gas cooling provided its particles interact with 3D phonons. The 
universal nature of 2D scattering at low energy of incoming particles makes 
these two assistance mechanisms to be robust to specific details of a particle- 
particle or impurity potential. The latter is only assumed to be short-ranged 
with the characteristic interaction radius ro being shorter than the DeBroigle 
wavelength. For example, an indirect exciton interacts with an impurity or 
another exciton via the electron- hole dipole moment: ed, directed along the 
normal to the bilayer. Hence, the radius of such a dipole interaction equals 
to the spacing between the electron and the hole layers: Tq ~ d. 

The Hamiltonian of particle phonon interaction can be written generally 



2 



as: 

Hx-ph = j j i^^{^W)5{z)Ti{p-p')u,{p') £p'd''fdz, (1) 

where p = {f,z), ip^ and ip are the particle creation and annihilation oper- 
ators and M is a crystal deformation induced by an acoustic phonon. Wave- 
lengths of relevant phonons are much larger than the width of the bilayer: d. 
The lattice deformation u can be expanded into the normal phonon modes 
as: 

""'^^ = ^ (2;^] (^^'^^(-^"^ + (2) 

where and bg are the phonon creation and annihilation operators of polar- 
ization s, p is the mass density of solid, Wsiq) = cq is the phonon frequency 
dispersion, which we assume to be isotropic and independent of phonon po- 
larization s. 

The exciton phonon interaction in GaAs crystal can be separated into 
piezoelectric and deformation potential parts. Lattice deformation in a piezo- 
electric crystal induces a polarization density: Pi = PijkdjUk 0, where Pijk is 
the piezoelectric tensor. This polarization interacts with the exciton dipole 
moment. In the limit qd <^ 1, the deformation potential for an exciton: B, is 
a sum of the deformation potentials for an electron and a hole taken at the 
same spatial point. The latter represents a change of the semiconductor gap 
due to the local compression caused by a phonon deformation. Combining 
the piezoelectric and deformation parts and expanding the crystal lattice de- 
formation in acoustic phonon modes we write the the exciton-phonon vertex 
in the Hamiltonian (|l|) as P: 

r.(g-)= (eg. + erfA,,^^^^)- (3) 

For a cubic GaAs crystal without the inversion center f3ijk = /3 if all i,j,k 
are different and zero otherwise. In the limit of large d the piezoelectric 
part dominates over the deformation potential but in the experiments ^ 
d ~ 50A and both exciton-phonon interaction terms are of the same order of 
magnitude: G ~ Anedp ~ lOeV^. 

Amplitude of phonon emission is given to the lowest order of the pertur- 
bation theory by the following matrix element: 

M!j{q) =< fqs\H,,\tO > (4) 
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between initial state of Bose gas 1^0 > with no phonons and the final state 
of Bose gas \fqs > with just one phonon specified by the momentum: q and 
the polarization s. We assume that the thermalization of the Bose gas due 
to particle-particle scattering is much faster than the slow cooling due to 
phonon emission. Thus, at any given time t the Bose gas is characterized 
by an effective temperature: T{t) . This temperature defines the total gas 
energy: E = E{T). The Fermi Golden Rule gives the probability of phonon 
emission per unit time and one needs to multiply it by the phonon energy: 
ujs{q) = cq, to find the total energy losses: 

^ = -X E cq\M:^m'm -e^- cq). (5) 

fqs 

Eq. (H) has to be averaged over Gibbs distribution of the initial states with the 
effective temperature T{t). Both the initial and the final states of the Bose 
gas are calculated in the interaction representation (see e.g.p]). Particles are 
confined to the 2D layer and the energy losses are proportional to the area 
of this layer. 

In the experiments [|I], 0] only a relatively low density of exciton gas has 
been achieved. Popov has shown that for a 2D dilute Bose gas there is a 
2D Berezinski-Kosterlitz-Thouless phase transition point: 

Tc = : r (6) 

gm log L 

that separates high-T almost ideal Bose gas phase from the low-T superfiuid 
phase. Actually, Popov theory is controlled by the large logarithm: 

L ^ - log{nrl) ^ log ^ (7) 

J- c 

where n is the particle density, Eq = h'^/r-QUi, g is the particle internal de- 
generacy [0]. For a Bose particle: g = 2S + 1, where 5" is the spin of particle. 
It was shown in Ref. p[ that electron and hole spins fiip rapidly due to the 
spin-orbit interaction. Thus, g = 4 ioi a. GaAs exciton. 

For 2D dilute non-ideal Bose gas one can distinguish three temperature 
regions. At high temperatures: T ^ TclogL, the ideal Bose gas is a good 
approximation. At intermediate temperatures: gTc/L -C T <^ Tc/ log L, an 
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Figure 1: Amplitudes of exciton-exciton scattering (left) and impurity scat- 
tering (right) accompanied by an emission of phonon - shown by a wavy 
line 



overwhelming amount of particles constitute a 2D BEqC with the density: 

= (l - ^) , (8) 

whereas a small fraction of thermal particles have a bare dispersion and Bose 
distribution with the chemical potential: fi ~ gTc/L |0. At low tempera- 
tures: T fi, a weak particle-particle interaction is crucial and the Bose 
gas becomes a Bose-liquid with quasiparticle excitations having Bogolubov 
sound like dispersion. The transfer of momentum to an impurity becomes 
here inefficient because sound-like quasiparticles ignore point-like impurities. 
Hot particles with a bare dispersion are rare due to an exponentially small 
Boltzmann factor. In the case of excitons in GaAs crystal the intermediate-T 
region hardly exists at all. 

Fig.l shows two scattering possibilities for a particle accompanied by an 
acoustic phonon emission. The left diagram shows a scattering on the second 
particle and the right diagram shows a scattering on impurity. First, we treat 
this problem in the high-T and intermediate-T regions where particles have 
the bare dispersion: e = /2m. 2D scattering amplitudes in both cases are 
isotropic and depend only on the total kinetic energy in the center of mass 
frame: £ (in the impurity case £^ = e) in the limit £ ^ 

F{n = -(^og-^) , 9 

m* \ £rQm J 

where m* = m/2 for particle-particle scattering and m* = m for the impu- 
rity scattering (see e.g. M). Both particle-particle interaction line and the 
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impurity line on Fig.l correspond to the scattering amplitude F{£). The 
total amplitude of assisted phonon emission is universal in both cases and is 
given by the following matrix element: 



M^(g-) = C iF{£)- + F{£ - cq)^\ 
\ cq —cq J 



efr.(,-)J^, (10) 



where C = 2 for the particle-particle scattering and C = 1 for the impurity 
scattering. We neglect the phonon momentum gj| transfer to the particle 
because in the low temperature limit: q\\ -C p, -C cm. 

Plugging Eqs.(p!0|,|3|) into Eq.(||) and taking the integral over the final state 
of the Bose gas we obtain the total energy losses per unit time. In the high-T 
region we get: 

dE 27lC^ f( ^ 1 rr^^u > ^JJ|d^q 



/ (— --) K(e,T)5(e-e'-cq)dede' 

J Woe En /e log En /e' / ^ ' ^ ^ 



dt hT^_^J \\og{Eo/e) \og{Eo/e')J '''' ^' (27r)3' 

(11) 

where K{e,T) = 2AnN{e){l + l/g) in the case of particle-particle scattering 
and K{e, T) = AnimpN{e)/2 in the case of impurity scattering. A is the total 
area of bilayer, riimp is the areal density of impurities and 

Nie) = \ (12) 

exp((e-/i)/T) - 1 

is the Bose-Einstein occupation number. It is also convenient to define a 
characteristic exciton-phonon temperature in the case of GaAs: 



x—ph 



^ -5°K. (13) 



\ 92 + {47idepy/15 



Combining the particle-particle and impurity contributions and using the 
ideal gas equation of state: E{T) = AnT, we find the total cooling rate: 

dT ,,,, , , , , 85 T3 

- = -(4(l + l/,). + n..,)--^. (14) 

Integrating the Eq.(0) we get the temperature relaxation law: T(t) ~ 1/v^. 
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In the intermediate-T region at g = 1 the coohng rate is enhanced by the 
stimulated scattering into the BEqC final states: 

In the low-T region the thermodynamic equation of state reads: E{T) = 

A(^{3)T^/7is'^, where s is the Bogolubov sound velocity: s = \J ji/m. In order 
to calculate the energy losses we apply a Bogolubov unitary transformation 
to the Hamiltonian (|l]): 

Hvh=-Yl 7rn{^''{P)^''{-P + miQ)Ui{q) + c.c. (16) 

pq ^^IPJ 

This Hamiltonian allows an emission of a phonon. The cooling rate in this 
case is also enhanced by the condensate stimulation: 

dt MTl^^X ({3)) ■ ^ ' 

For experimental realization of an exciton cooling the most relevant is 
Eq.(|T^). Integrating it and inserting L ^ 6 and n ^ Uimp we find the overall 
cooling time tc required to reach the phase transition point. Here one can 
distinguish two cases: i) cooling of exciton gas with constant density e.g. 
sustained by photoexcitation and ii) cooling of decaying exciton gas with 
n{t) = rig exp(— t/r^), where is the exciton recombination time. We find: 

= (18) 

where Tc(n) is the BKT temperature as a function of the exciton density n 
(^, C ~ 10 in the case i) and C ~ 30 in the case ii). Note that C does not 
depend on exciton mass m and in the case ii) the best condition for reaching 
Tc occurs at t = r^/S. Eq.([T^) defines a frontier in the bilayer parameter 
space: (n, r^) or equivalently (Tc{n),Tr), separating the two kinetic phases - 
one that can condensate and the second that remains above Tc during the 
exciton life time r^. Fig. 2 shows this borderline for the case of a thin GaAs 
bilayer. 

This work was supported by RFFI and INTAS. 
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Figure 2: Kinetic phase diagram showing the posibihty to reach a Bose 
Einstein quasi-condensate phase in the GaAs indirect exciton bilayer. 
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